Observable to explore high density behaviour of symmetry energy 



O 
(N 

D 
00 



> 

m 
o 



Aman D. SoodlB 

^ SUBATECH, Laboratoire de Physique Subatomique et des Technologies Associees 
University of Nantes - IN2P3/CNRS - Ecole des Mines de 
Nantes 4 rue Alfred Kastler, F-44 072 Nantes, Cedex 03, France 
(Dated: September 29, 2011) 



Introduction 

The nuclear matter equation of state (EOS) 
of asymmetric nuclear matter has attracted 
a lot of attention recently The EOS of 
asymmetric nuclear matter can be described 
approximately by 

E{p,5) = Eo{p,5 = Q)+EsyM5'' (1) 



where 5 



Pn~Pp 



^ is isospin asymmetry, Eo(/c, 
5) is the energy of pure symmetric nuclear 
matter, and Esym(/3) is the symmetry energy. 
The symmetry energy is important not only 
to the nuclear physics community as it sheds 
light on the structure of radioactive nuclei, 
reaction dynamics induced by rare isotopes, 
but also to astrophysicists since it acts as a 
probe for understanding the evolution of mas- 
sive stars and the supernova explosions. Al- 
though Esym at saturation density is known 
to be around 30 MeV its value at higher densi- 
ties is poorly known. Experimentally, symme- 
try energy is not a directly measurable quan- 
tity and has to be extracted from observables 
which are related to symmetry energy. There- 
fore, a crucial task is to find such observables 
which can shed light on symmetry energy. In 
this paper, we aim to see the sensitivity of col- 
lective transverse in-plane flow to symmetry 
energy at low as well as high densities and also 
to see the effect of different density dependen- 
cies of symmetry energy on the same. The var- 
ious forms of symmetry energy used in present 
study are: Esym k (""), Egym fx (w)°'^, and 
Esym oc where u = The different 



density dependencies of symmetry energy are 
shown in Fig. 1. The present study is car- 
ried out using IQMD model which is described 
briefly in the next section. 



The model 

In IQMD ll model the hadrons propagate 
using Hamilton equations of motion: 



d'rl d(H) dp. 



dpt 



dt 



d{H) 
dfi 



where (2) 



{H) = {T) + {V) 

i i j>i 

xfjif',p',t)dT'df'dpdp'. (3) 

The baryon potential V^^ , in the above rela- 
tion, reads as 



Skyrrne ~ Yukawa ~ Coul 

(4) 
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Where the potential terms in eq. 4 represent, 
respectively, Skyrme, Yukawa, Coulomb, mo- 
mentum dependent interaction, and symme- 
try potential. 



Results and Discussion 

There are several methods used in the liter- 
ature to define the nuclear transverse in-plane 
flow. In most of the studies, one uses {px/A) 
plots where one plots {px/A) as a function of 
^c.m./^beam' Using a linear fit to the slope. 
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FIG. 1: Different forms of symmetry energy. 



one can define the so-called reduced flow (F). 
Alternatively, one can also use a more inte- 
grated quantity "directed transverse in-plane 
flow (p^^'^)" which is defined as Q: 



(5) 



where Y{i) and Px(i) are the rapidity distribu- 
tion and transverse momentuni of the ith par- 
ticle. In this definition, all rapidity bins are 
taken into account. It therefore presents an 
easier way of measuring the in-plane flow than 
complicated functions such as (p^/A) plots. 

In Fig. 2 we display < ^ > as a function of 
Yc.m./y|3cam final time (left panels) and 
the time evolution of < p'^^^ > (right pan- 
els) calculated at 100 (top panel), 400 (mid- 
dle) and 800 MeV/nucleon (bottom) for differ- 
ent density dependencies of symmetry energy. 
Solid, dash dotted, and dotted lines represent 
the symmetry energy proportional to p, 
and p^, whereas dashed line represents calcu- 
lations without symmetry energy. Comparing 
the left and right panels in Fig. 2, we find 
that both the methods show similar behavior 
to symmetry energy. For example, at incident 
energy of 100 MeV/nucleon for Esym oc p"-^, 

< Px^^ > = 0. Similarly, the slope of < ^ > 
at midrapidity is zero. We also find that 
the transverse momentum is sensitive to sym- 
metry energy and to its density dependences 
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FIG. 2: Left panel: < ^ > as a function of 

Yc.m./yjjgaui and < p^^^ > at different incident 
energies. Panels and lines are explained in the 
text. 

Fi(u), F2(it) and F3(u) in the low energy re- 
gion (100 MeV/nucleon) only. At energies 
above Fermi energy, both the methods show 
insensitivity to the different symmetry ener- 
gies. This could be because (i) the repulsive 
nn scattering dominates the mean field at high 
energies, (ii) As explained in ref Q that the 
effect of Esym on nucleons during the initial 
stages of the reaction (0-30 fm/c) decide the 

fate of final state < p^^"^ >. At low energies 
duration of high dense phase will prevail for 
a longer duration alllowing the E^ym to affect 
the flow. Whereas at high energies, although 
maximum density reached will be larger but 
the duration of high density phase will be very 
small not providing Egy^ with sufficient time 
to make an impact on flow. 
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